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Time-dependent density functional theory, proposed recently in the context of atomic diffusion 
and non-equilibrium processes in solids, is tested against Monte Carlo simulation. In order to assess 
the basic approximation of that theory, the representation of non-equilibrium states by a local 
equilibrium distribution function, we focus on one-dimensional lattice models, where all equilibrium 
properties can be worked exactly from the known free energy as a functional of the density. This 
functional determines the thermodynamic driving forces away from equilibrium. In our studies of 
the interfacial kinetics of atomic hopping and spin relaxation, we find excellent agreement with 
simulations, suggesting that the method is useful also for treating more complex problems. 
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I. INTRODUCTION 

Establishing a link between macroscopic laws of dif- 
fusion and relaxation with a microscopic master equa- 
tion for atomic degrees of freedom has remained a funda- 
mental problem in non-equilibrium statistical mechanics. 
Moreover, numerous examples exist in metallurgy, semi- 
conductor device technology, glass science and polymer 
science, where control over the time-development of mi- 
crostructures is crucial in the design of materials with 
special mechanical, electrical and magnetic properties. 
Hence there is a need also from a practical viewpoint to 
derive tractable kinetic equations incorporating specific 
materials properties, so that processes like nucleation, 
spinodal decomposition and magnetic relaxation can be 
treated in a realistic manner. 

The simplest approach is to study mean-field kinetic 
equations, derived from the master equation by neglect- 
ing any atomic correlation effects pi El 9]. Such equa- 
tions suffer from the obvious drawback that their station- 
ary solutions yield the mean-field phase diagram, which 
may differ even qualitatively from the correct phase di- 
agram. Several routes for improvement have been pro- 
posed in the literature, including the path probability 
method (4|, || , effective Hamiltonian methods B or time- 
dependent density functional theory (TDFT) J?], |, |). 
The latter approach implements the idea of local equilib- 
rium and leads to thermodynamic driving forces which in 
principle are derived from an exact free energy functional. 
Density functional theories are normally formulated for 
continuous fluid systems fl(i[| , but adaptation to discrete 
lattice systems is straightforward |ll], |l2). In this way 
one can obtain generalized mean-field kinetic equations 
for single-particle or single-spin densities which in princi- 
ple are consistent with the exact equilibrium properties. 



If necessary, additional approximations with respect to 
equilibrium quantities can be carried out in a separate 
step. 

The above-mentioned kinetic theories mostly focus on 
purely dissipative processes in discrete lattice systems. 
For fluid systems, the derivation of nonlinear transport 
equations for hydrodynamic variables is a more com- 
plicated subject. Some generalizations of Moris' well- 
known projection operator technique and mode coupling 
approximations to situations far away from equilibrium 
Jl3| . [lJJ have been applied, for example, to problems of 
nonlinear hydrodynamics, the glass transition ]l5[ | and to 
granular flows fLq| . 

Our aim in this paper is to apply the TDFT scheme 
to purely dissipative "conserved" atomic migration 
and "non-conserved" spin dynamics processes in one- 
dimensional lattice models for which the exact free en- 
ergy functional is known. This enables us to separate 
out and to test the local equilibrium assumption against 
Monte Carlo simulation. Specifically, we study the tem- 
poral evolution of domains with different ordering, and of 
the associated interface. It is demonstrated that in these 
problems the TDFT shows remarkable quantitative ac- 
curacy, suggesting that this method may be useful also 
under more general conditions. 

In section 2 we briefly recall our basic approach. Sec- 
tion 3 starts out with an exact free energy functional for 
a one-dimensional lattice and provides expressions for lo- 
cal correlators in terms of particle densities. With these 
results we arrive at a closed system of kinetic equations 
on the single-particle level. Following the classification 
by Hohenberg and Halperin Jl7|] these equations take the 
form of generalized "model B" or "model A" equations in 
cases of a conserved or a non-conserved order parameter, 
respectively. Solving these equations, we subsequently 
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discuss the time-evolution of an initially sharp interface 
between two differently ordered domains in both of these 
cases (section 4) . Excellent agreement with Monte Carlo 
simulations is found, in contrast to ordinary mean-field 
(MF) theory which produces substantial deviations. 



II. OVERVIEW OF TIME-DEPENDENT 
DENSITY FUNCTIONAL THEORY (TDFT) FOR 
STOCHASTIC LATTICE SYSTEMS 

A. Atomic hopping 

Let us begin with hopping of particles on a lattice 
of equivalent sites i, which are either simply occupied 
(rii = 1) or vacant (n, = 0), so that occupation num- 
bers satisfy nf = rii. The hopping process is described 
by a master equation for probabilities -P(n, t) of finding 
an occupational configuration n = {rii} at time t. As 
elementary steps we assume moves of a single particle 
from an occupied site to a vacant nearest neighbor site. 
The associated rates Wi tk (n) for adjacent sites i and k 
to exchange their occupation satisfy the detailed balance 
condition with respect to a given lattice gas Hamiltonian 
H(n). 

A detailed description of TDFT is found in Ref. |§. 
Hence we need to recall only the main steps, and add 
some remarks as to their physical content. The basic 
approximation is to replace the distribution P(n, t) by 
the local equilibrium distribution 
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where deviations from equilibrium are represented in 
terms of time-dependent single-particle fields h(<) = 
{hi(t)}. Z(t) is a normalization factor, which at equilib- 
rium (h(i) = 0) reduces to the canonical partition func- 
tion. Requiring self-consistency on the single-particle 
level allows us to eliminate h(i) in favor of mean occu- 
pation numbers p(t) = {pi(t)} with pi(t) — (ni) t , where 
(. . . ) t denotes an average with respect to the distribution 
(|l|). In this way a closed system of equations for p(t) can 
be derived. 

To carry through this program we start from the equa- 
tion of continuity which follows directly from the original 
master equation. Replacement of exact averages by local 
equilibrium averages gives 



dpi(t) 
dt 



kit 



0, 



(2) 



with known expressions |q] for the current ji tk (n) from 
site i to site k in terms of Wi }k (n). Notice that at any 
instant of time the exponent in ([!]) describes an inho- 
mogeneous lattice gas which involves a spatially varying 
single-particle potential h(t). Hence, calculation of av- 
erages from (Q) is precisely the kind of problem treated 



by density functional theory (DFT) in classical statisti- 
cal mechanics. There, one considers a class of systems 
with fixed interactions and arbitrary single-particle po- 
tentials, specified here by H(n) and h, respectively. Av- 
eraged occupation numbers and correlation functions are 
determined from derivatives of a free energy functional 
F(p) associated with the Hamiltonian H(n). Specifically, 
p(t) is determined by the set of equations 



hi{t) + m(p(t)) = iHoU 



(3) 



with ptot the overall chemical potential, and 

W (p) = dF(p)/d Pl , (4) 

the local chemical potential as functional of p. Much ex- 
perience has been gained during the last two decades how 
to construct F(p) from a given Hamiltonian H(n). In the 
subsequent considerations we therefore assume F(p) to 
be known. Since in the framework of DFT occupational 
correlation functions are functionals of p, we can formally 
regard ([^) as the desired closed set of equations for p(i). 

In order to make this procedure explicit and to estab- 
lish a connection with thermodynamic driving forces, we 
again recall Ref. || where it is shown that the average 
current can be written as 



{ji,k)t = M i<k {t)[Ai{t) - A k {t)\. 



The quantities 



Mt) = exp[/3w (p(t))] 



(5) 



(6) 



are local activities, whose discrete gradient (along the 
bond connecting i and k) plays the role of a thermody- 
namic force that drives the current. The quantity 

M hk (t) = - (w itk (n) exp [{3{h t (t)n t + h k {t)n k )]) t , (7) 

where Mi jk (t) = M k ,i(t), is a mobility coefficient that 
depends on the actual nonequilibrium state. Further dis- 
cussion of (0) simplifies when we choose the hopping rates 
Wi ik (n) such that they depend only on the energy in the 
initial state, i. e. 



iu»,fc(n) = a[m(l - n k )e^ Hl + n k (l ~ rii)e 



JiH k 



(8) 



The first term describes hopping from i to k, with a ther- 
mally activated rate determined by the interaction energy 
Hi of a particle at site i with its environment, a is some 
bare rate constant. The reverse hopping process is de- 
scribed by the second term in (|§|). With this expression 
for 10^(11), one can show || that Eq. (0) transforms into 



M iik (t) = a((l- ni )(l-n k )) t . 



(9) 



At this stage, Mi jk (t) does no longer explicitly depend 
on h. Physically, equation (^) tells us that the mobility 
coefficient based on (||) is given by the nearest-neighbor 
vacancy correlator. 
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It should be kept in mind that in this TDFT -scheme 
all deviations from equilibrium are described in a mean- 
field manner in the sense that the underlying distribu- 
tion function (Q) deviates from the canonical distribution 
merely by single particle terms. Relationships between 
occupational correlators and densities p(t), which enter 
this theory, are local in time and are given by the equilib- 
rium theory. This implies the assumption that correla- 
tors relax fast to their local equilibrium values, compared 
with time scales characterizing the evolution of p(t). Or- 
dinary kinetic mean-field theory is recovered when we 
use mean-field expressions for /Xj(p) and replace (^|) by 
M ifk F i t ) = a(l-p*(<))(l-Pfc(£)). By contrast, in TDFT 
the local chemical potential appearing in (g) is defined 
by the exact chemical potential functional so that (|^) 
together with (^) describes relaxation towards the exact 
equilibrium state. Moreover, the expression (|J) for the 
mobility preserves local correlation effects in the jump 
dynamics. 



B. Spin relaxation 

To exemplify the dynamics of a non-conserved order 
parameter, we study spin-relaxation in a kinetic Ising 
model. Elementary transitions in the underlying mas- 
ter equation are supposed to be individual spin flips 
<ji — * — (7j, where er^ = ±1. By Wi(cr) we denote the asso- 
ciated rate in an initial spin configuration er. The local 
equilibrium distribution P^ oc ^(<r,i) is analogous to (Q). 
It involves the Ising Hamiltonian H(cr) supplemented by 
time-dependent magnetic fields h(t), which couple to the 
spins in the form — hi(t)<Ji. As shown in the Ap- 
pendix, the equations of motion read 



dt V o{ai) t 



with kinetic coefficients 



T i {t) = 2{w i {<j)e-f }h ^) t 



(10) 



(11) 



F is the intrinsic free energy functional associated with 
the exchange interaction, and h an overall constant mag- 
netic field. Equation ( |Io| ) again displays the exact ther- 
modynamic driving force in the spirit of TDFT. It can 
be regarded as a generalized "model A" equation in the 
classification by Hohenberg and Halperin [ fl7| , whereas 
Eqs. (||), (||) and (||) constitute generalized "model B" 
equations. Note that sufficiently close to equilibrium one 
can ignore h(t) in (jl^) and linearize the sinh-term in ( |To| ) 
to obtain d{ai) t /dt ~ -2l3{w l {rr)) eq {dF/d{a l ) t -h). The 
kinetic coefficient is then simply given by the equilibrium 
spin flip rate {wi{cr)) eq . 



C. Consistency with thermodynamics 



The total free energy decreases monotonously with time 
until the equilibrium condition is satisfied. For hopping, 
the rate of change of the free energy is given by 

dF ^ dF d Pi 
~dt~2-j ~dp~~dt 

=-k B Tj2Mi, k Xi(e Xi - e x <°) 



i . k 



knT 



Y,Mi, k { Xi - x k ){e Xi - e Xk ) < 0, (12) 



i.k 



where we have used (0), (|J) and (^) together with the 
abbreviation d((3F/dpi) — Xi, and Mi >k = M k ^. Cur- 
rents through the boundaries of the system are supposed 
to be zero. The inequality in ( |l2| ) arises from M,^ > 0, 
see (0), and from (x — y){e x — e y ) > for x ^ y. The 
equality sign in ( |l2|) holds if and only if Xi = x k for all i 
and k, which means that jii — const. 

Similarly, for the kinetic Ising spin model, the total 
free energy including the coupling to the external field h 
satisfies 



d_ 

dt 



^-A^ffij - -k B T r fe^fe si nh x k < 0, (13) 



where x k — j3(dF/d(cr k )t — h). The inequality follows 
because > (sec (O)) and zsinhx > for x =/= 0. 
Equation ( |l3| ) becomes an equality if dF/d(<Tk)t — h for 
all k. 



III. ONE DIMENSION: EXACT FUNCTION ALS 

To test the local equilibrium distribution ([!]) it is de- 
sirable to avoid any approximation with respect to static 
properties. This can be achieved by using exact free 
energy functionals, which arc available for certain one- 
dimensional systems [[l8[ |l9[ p0[ . 



A. Atomic hopping 

For a lattice gas with nearest neighbor interactions on 
a linear chain of sites i; 1 < i < M; with occupied bound- 
ary sites at i = and i = M+l, the free energy functional 
reads (lil 



M M-l 

i=0 i=0 



- (1 - Pi) ln(l- 



.n=l 

(14) 



Finally it is easy to show that our evolution equations 
are consistent with the second law of thermodynamics: 



(n) 

where V denotes the interaction constant, and p i+1 i with 
n = 1, ...4 are the two-point correlators for the four 
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possibilities of particles and holes on site i and site i + 1, 



B. Spin relaxation 



Pi+i,i = ((! - n i+1 )ni) =Pi-p, 



(i) 



(4) 
°i+l,i~ 



= ((1 - n»+i)(l - ni)) 



Pi+la' 
l- P . 



i (1) 
Pi+1 + «+!,, 



From the techniques of Ref. p9[ it follows that 



(1) (4) 
Pi+l,i Pi+l,i 



P 



(2) _(3) 
i+lA Pi+XA 



-pv 



(15) 



(16) 



a relation, which is equivalent to the quasi-chemical 
approach. Given these relations, it turns out that 
dF/dp^n = 0. This suggests that F may be minimized 
also with respect to correlators in cases where these can- 
not be calculated explicitly. 

The boundary conditions for the correlators are p^Q = 
Pi and j4f+i m = Pm- F° r 1 < i < M — 1, combination 
of ( |l5| ) with ( |i~6| ) yields a quadratic equation for f . 
In this way the representation of ^ as functionals of 

p is completed. The fact that j only depends on Pi+i 
and pi clearly is a special feature in one dimension. 

The kinetic equations derived in section 2 are now com- 
bined with (0). Evidently, from (§), 

M M+1 (t) = apjjy, (17) 
while the local chemical potential is found to satisfy 

Pi 



(2) (3) 
, Pi+l,iPi,i-l , 

m ,, s ' — — In ■ 



„(4) (4) 
Pi+l,iPi,i-l 



1 



(18) 



From (§), @ and @ we obtain for the current 

(ji,i+l)t = 



p?L (i 



Pi J (2) 
Pi+1, 



* Pita 

Pjjkj+l (1 -Pi+l) (3) 

(4) Pj+l,j 

Pi+1 n y ' 



(19) 



In these last equations (|17[)-(|19|), densities and correla- 
tors are local equilibrium quantities. The final form of 
our kinetic equations as a nonlinear set of differential 
equations for pi (t) emerges when we reexpress p^\ i in 
terms of Pi+i and pt in the way described above. 

For comparison we also consider the ordinary mean- 
field equations. These are obtained by factorizing all cor- 
relators in ( |l~8| ) and (|l9|), for example p\V lt — Pi+ipt- 
The mean-field current is then found as 



Uifi+i)t = " {Pi - Pi+i + K\Pi-xPi{L - Pi+i) 

-(1 - Pi)pi+lPi+2}} 

with K = exp(fiV) - 1. 



(20) 



Next we turn to spin relaxation in the linear Ising 
model. Rather than using (10) we immediately choose 
transition rates 



a / 7 

u'i{cr) = — (l - -fJ^CTj+i + <7j_i) ) (1 - 5a. l 



(21) 



and start from the evolution equations for single spins as 
given in the original work by Glauber EM, 



dt 



o-i-i)t 



+ Y« Cr i+l (J i>t + (<W-l)t) 



(22) 



Here, 7 = tanh2/3J and 6 = tanh/3/i, where J > 
denotes the (ferromagnetic) nearest-neighbor exchange 
coupling and h a constant external magnetic field. It 
is well known that for h = these equations become 
linear and easily soluble. By contrast, for h 0, the ap- 
pearance of correlators {(Ti+i(Ti)t in ( [22] ) prevents us from 
obtaining an exact solution. Using the well-known repre- 
sentation of Ising spin variables by occupation numbers, 
<7j = 2ni — 1, and vice versa, we can treat the correlators 
(ai+iOj) in perfect analogy to (ni+irii}. In particular, 
Eq. ( jlq ) with J = 4U transforms into a quadratic equa- 
tion for (<7i±i<Ji), whose solution, expressed in terms of 
(cy-i} and (o-j), is substituted into (p2|). This yields our 
TDFT-equation of motion for spins. Likewise, we ob- 
tain from ( |l6| ) the free energy as functional of the spin 
density, which could be used in (jl0|). 



IV. APPLICATION TO INTERFACIAL 
KINETICS 

We now apply the TDFT to problems of the time- 
evolution of an initially sharp interface between differ- 
ently ordered domains on a linear chain. Our purpose 
is to present a quantitative comparison with both Monte 
Carlo simulation and simple MF-theory with respect to 
density profiles, spin-density profiles and the respective 
correlators. 



A. Atomic hopping 

The length of the chain is taken as M = 10 3 . As 
mentioned before, boundary sites have fixed occupation 
Pa = Pm+i = 1- Symmetrical initial conditions at 
t = are chosen such that we have a vacant region cen- 
tered around the midpoint of the system, Pi(0) = for 
250 < i < 750, and complete occupation in the com- 
plementary space. For t > 0, the initially sharp density 
profile will progressively broaden due to diffusion. This 
is shown in Fig. [l] for < i < 500 in the case of a re- 
pulsive interaction with [3V — 3. Generally, the shape 
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FIG. 1: Comparison of time-dependent density profiles pi 
in the case of hopping dynamics with repulsive interaction 
13V = 3, obtained from different methods, a) TDFT (lines) 
and MC simulation (data points) b) Kinetic MF-theory. 



of profiles depends on how the interaction enters the el- 
ementary hopping rates. Our choice (ph implies that in 
regions with densities p > 0.5 a particle next to a vacant 
target site has a large chance to be repelled by another 
particle and hence will assume a large jump rate. By con- 
trast, the repulsion will be less active in dilute regions. 
This explains the asymmetry of the profiles in Fig. 0a, 
with a steep drop towards the empty region. The main 
conclusion from Fig. |l|a is the perfect agreement between 
profiles from TDFT, shown by the full lines, and from 
Monte Carlo (MC) simulation (data points) @. To get 
smooth profiles from simulation, we took averages over 
10 4 Monte Carlo runs. By contrast, the MF-profiles in 
Fig. |l|b are more symmetric and deviate significantly from 
those in Fig. |l|a. 

For diffusion processes on (continuous) length scales x 
and time scales t much larger than the elementary hop- 
ping distance and residence time, we expect the density 
to depend only on the scaling variable r\ — xj pro- 
vided the initial conditions can be expressed in terms of rj. 
This is verified in Fig. |^a which shows master-curves p(rj) 
obtained from the profiles in Fig. [I] for different times. In 
this analysis the origin of the x-axis is chosen to coin- 



FIG. 2: a) Density profiles shown in Fig. 1 for different times 
against the scaling variable rj — x/(2yi). The length x in 
units of the lattice constant and t in units of At, see foot- 
note pi| . b) Concentration-dependent diffusion coefficients 
D(p) extracted from the master-curves of Fig. 2a by the 
Boltzmann-Matano method. (The normalization factor Do 
is the single-particle diffusion constant for infinite dilution.) 

cide with the initial density drop at i = 250. As ex- 
pected from Fig. |l|, the TDFT master-curve, in contrast 
to MF-theory, practically coincides with the Monte Carlo 
master-curve. 

These results can be analyzed further by the 
Boltzmann-Matano method |25| , which assumes a diffu- 
sion equation of the form dp/dx — d/dx(D(p(x))dp/dx) 
to hold. From the master-curve p(r]) the concentration- 
dependent diffusion coefficient D{p) can be deduced ac- 
cording to 

where rj(p) is the inverse function of p{rf). The integral 
can be calculated accurately from the profile of Fig. ||a up 
to p ~ 0.9, and the results for D(p) are shown in Fig. ||b. 
Using the MF profile, we recover the p-dependence of the 
mean- field diffusion constant. This quantity is calculated 
easily by separating from the current (|2^) a factor Pi+i — 
Pi, i. e. a discrete gradient of the density, and identifying 
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500 



FIG. 3: Correlators (nj+iru)t at t — 10 4 from MC-simulation 
compared to correlators computed from TDFT and MF- 
theory, using the same Monte Carlo density profile pi (upper 
curve) as input. The TDFT-correlators are indistinguishable 
from MC-correlators in this plot. 



the result with Ficks' law. One obtains 



D MF {p)=D Q {l + K[p 2 + Ap(\-p)]), (24) 



-l 



V 1 






x X \ a ) 








\ ' t 1 

\ 10 \ 

\ioo \ \ 
Mm \ u 




\ 700 


400 






1000\ ' ' » . 







10 20 30 40 50 



-1 







. V\\ 


b) 






\ io 

20 i 

40 


t= 1 










-.100 

:ooo ___ 









10 20 30 40 50 



with Do = aa 2 , a being the lattice spacing. The expres- 
sion ( p4[ ) shows a broad maximum around p — 2/3, which 
reflects the average effect of the repulsion of particles. 
The TDFT-diffusion constant, however, shows a much 
sharper maximum. Moreover, when p becomes small, it 
approaches the value Dq more rapidly, and thus gives 
rise to the steepening of the density profile in the regime 
p < 0.4, as observed in Fig. ^a. This is a correlation ef- 
fect: In a dilute system, a fast hop of a particle due to the 
repulsion by a neighboring particle is a rare event because 
nearest neighbor pairs get suppressed, (ni + in,i) < pi + \Pi, 
and hence diffusion is slowed down relative to the MF- 
prediction. This argument is supplementary to our pre- 
vious discussion of Fig. [j]a. 

At this point we remark that such correlations induced 
by the repulsion of particles are taken into account to 
a certain extent even by MF-theory when applied to a 
two-sublattice structure. Density profiles and effective 
diffusion coefficients calculated in this way in a previ- 
ous study Q indeed are similar to those of the present 
TDFT calculation shown in Fig. ||a and ||b, respectively. 

Because of the important role played by the correla- 
tors in TDFT it is of interest to make a direct compari- 
son with correlators obtained from Monte Carlo simula- 
tion. Fig. [| exemplifies perfect agreement between those 
of TDFT and simulation, whereas MF-correlators, cal- 
culated here as product pi + \pi of the simulated densities, 
are significantly larger when the densities are small. 



FIG. 4: Comparison of time-dependent spin density profiles 
in the Glauber model with /3J = 2 and /3/i = —0.1, obtained 
from different methods, a) TDFT (lines) and MC-simulation 
(data points) b) Kinetic MF-theory 

B. Spin relaxation 

In our study of "non-conserved" dynamics in the one- 
dimensional Glauber model |2l[] we choose a chain of 
length M = 10 2 and fixed upward spins at the bound- 
aries, o"o = cm+i = I- Our initial condition at t = now 
is Ui = — 1 for 40 < i < 60 and oi = +1 for the remaining 
spins. Notice that in the case h = simple MF-theory in 
the spirit of this work becomes exact because the correla- 
tors in (|2^) drop out. To depart from this trivial situation 
we introduce a small field with (3h = —0.1 which favors 
downward spin orientation. Spin density profiles in the 
region < i < 50 for (3 J — 2 at different times t > 
are presented in Fig. ^a, where the full lines correspond 
to TDFT, and data points to simulation. The agreement 
is very good, although not perfect. Generally, the spins 
in the interior of the system relax towards the equilib- 
rium in the external field. Spins near the boundary are 
expected to relax towards an equilibrium profile which 
decays from the boundary (do = 1) towards the interior 
(<7i ~ —1) on a length given by the correlation length 
£. For the temperature considered, £ ~ 5 a. During the 
course of this relaxation, simple MF-theory, based on a 
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Under the ultimate aim to develop theoretical tools for 
a description of phase transformation processes in real 
materials, several extensions of the present work are nec- 
essary. First of all, one needs reliable approximations for 
the free energy functional in higher dimensions. For two- 
dimensional lattice systems, a step in this direction has 
been taken recently p4j , which was based on an exten- 
sion of the techniques in Ref. ]l9| . Secondly, one would 
like to treat multicomponent systems. In that case, local 
equilibrium distributions of the type (Q) may be insuffi- 
cient to describe interdiffusion currents related to non- 
diagonal Onsager coefficients ||. To incorporate such 
effects into the TDFT-scheme is an open question which 
deserves further study. 



FIG. 5: Correlators (<7i+i<7,:)t at t = 10 2 from MC-simulation 
compared to correlators computed from TDFT and MF- 
theory, using the same MC spin density profile {<n)t as in- 
put. The TDFT-correlators are nearly indistinguishable from 
MC-correlators in this plot. 



factorization of the last two terms in (|22|), gives quite 
different results (Fig. First of all, the overall process 
is much faster than in Fig. |]a. Second, within the region 
of the initial upspin domain it predicts a constantly de- 
creasing plateau which is not observed in Fig. [|a. The 
origin of these failures of MF-theory becomes clear when 
we look at Fig. ||: Monte Carlo- and the almost identical 
TDFT-correlators (er i+ iCTi)t stay close to unity through- 
out the system, in contrast to the MF-factorization, and 
stabilize the respective spin configuration. Hence the re- 
laxation process progresses only by a successive broad- 
ening of the interfacial region between the upspin and 
downspin domains and not by a decaying plateau. 



V. SUMMARY AND OUTLOOK 

Applying a local equilibrium approximation to the 
master equation for atomic or spin configurations, ki- 
netic equations for particle or spin densities were derived, 
which are compatible with the exact thermodynamics. 
The derivation was largely based on concepts from den- 
sity functional theory. Kinetic equations obtained have 
the form of generalized "model B" or "model A" equa- 
tions in the language of Ref. |17[ ], where thermody- 
namic driving forces originate from the exact free energy 
functional. This "time-dependent density functional" 
(TDFT)-scheme is tested against Monte Carlo simula- 
tions for both a one-dimensional hopping model and the 
Glauber model, where the exact free-energy functional is 
known. Studying the dynamics of the interface between 
different domains, the TDFT yields excellent agreement 
with simulations with respect to density or spin density 
profiles and local correlation functions. The success of 
this theory appears to be a consequence of the fast relax- 
ation of correlators towards their local equilibrium values. 
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Appendix: Derivation of generalized "model 
A" —equations 



The derivation of Eq. ( 10 ) proceeds in steps with some 



similarity to Ref. |8j]. In the present "non-conserved" 
case we start from the master equation for single spin 
flips, see e.g. Ref. pl| , to obtain the time derivative 
of single-spin averages. Exact averages are in turn ap- 
proximated by averages (■■■)* with respect to the local 
equilibrium distribution p(' oc )(er, t), which has the same 
form as ([l]) apart from a sign change in the second term in 
the exponent. (This is because the auxiliary fields hi(t) 
in (|l|) have the meaning of effective site energies, while 
they are taken here as effective magnetic fields.) In this 
way we arrive at 



d((Ji 



dt 



= -2(wi(cr)ai 



(25) 



The summation over all <x in the definition of the average 
on the right-hand side of (^5|) involves a summation over 
(Tj = ±1, which we treat with the help of the detailed 
balance condition: 

2J exp[-f3(H((r) - hi(Ti)]wi(cr)ai 

= ^e- /5H(CT) w 4 (cr)sinh/3/i 4 . (26) 



In the last step we have used sinh ((3hi<Ji) = <Ji sinh(f3hi) 
and of — 1. To restore the expression for P( loc \cr,t) 
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we multiply and divide (Eq) by exp(/3hi). Finally, it fol- 
lows from the form of P^°^(<r,i) that single-spin aver- 
ages and the fields h(t) are connected by hi(t) + h = 
dF/d(<Ji)t, which is analogous to (||). Here, F is the 



intrinsic free energy as a functional of the spin density. 
Combination of these results with ( p5|) yields ( |l0|) and 
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